Quantitative analysis of the debonding structure of soft adhesives 
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We experimentally investigate the growth dynamics of cavities nucleating during the first stages of 
debonding of three different model adhesives. The material properties of these adhesives range from 
a more liquid-like material to a soft viscoelastic solid and are carefully characterized by small strain 
oscillatory shear rheology as well as large strain uniaxial extension. The debonding experiments are 
performed on a probe tack set-up. Using high contrast images of the debonding process and precise 
image analysis tools we quantify the total projected area of the cavities, the average cavity shape 
and growth rate and link these observations to the material properties. These measurements are 
then used to access corrected effective stress and strain curves that can be directly compared to the 
results from the uniaxial extension. 
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I. INTRODUCTION 

When soft adhesives are detached from rigid surfaces, 
the incomprcssibility of the material combined to its ex- 
treme deformability leads to complex deformation pat- 
terns involving the formation of air fingers and cavities 
[D-Q. The details of these patterns depend markedly on 
the material properties and often evolve towards a fibril- 
lar structure of highly stretched polymers which eventu- 
ally fail by fracture or detach from the surface @. The 
criteria for the onset of the initial elastic or viscous insta- 
bilities have been known for some time @, H[ and several 
experimental studies have focused on fingeri ng instabili- 
ties [sj-TTojl . on the cavitation criteria lllM l5l|. cavity nu- 
cleation rate 0, [l7| or growth rate [18| . However the 
transition from growth of individual cavities to the collec- 
tive growth of a population of cavities under the applied 
stress, leading to elongated walls between cavities, also 
called "fibrils" and to eventual detachment, has received 
much less attention experimentally 0, [l?} . Some theo- 
retical papers have been published on collective growth 

mill. 

Up to date, it remains difficult to relate the observed 
patterns to the rheological properties of the soft adhe- 
sives, mainly due to the lack of precise experimental 
characterization of the 3D structures and of the mate- 
rial deformation during the debonding process. Because 
the processes are dynamic, powerful 3D scanning tech- 
niques, such as confocal microcopy, are too slow and one 
has to rely on classical 2D imaging limited by its depth 
of field. Proper identification of the cavity borders in an 
automatic and reliable way it is not a trivial task and re- 
quires good quality well contrasted images and adapted 
imaging software tools. Yet, this information, albeit sta- 
tistical in nature, is essential if one wishes to gain more 
insights on the debonding process and be able to compare 
experiments with results from numerical simulations. It 
is also a necessary ingredient to understand which rheo- 



logical properties of the material determine the debond- 
ing patterns and, eventually, the adhesion performance, 
the important parameter for practical applications. 

In this paper we have performed careful experiments 
yielding high contrast images of the cavities nucleated 
in the early stages of debonding during a probe tack 
test. We have developed precise image analysis tools to 
characterize quantitatively and in a statistically signif- 
icant way the size, shape and overall projected surface 
of the cavities. Using model materials with well known 
rheological properties spanning from viscoelastic liquids 
to viscoelastic solids, we will present detailed measure- 
ments of the growth dynamics of cavities, including the 
total projected area, the average cavity shape and their 
growth rate. These measurements give access to a cor- 
rected true stress and strain which can then be quanti- 
tatively compared with material properties in shear and 
uniaxial extension. 



II. MATERIALS AND METHODS 

A. Sample preparation 

Model acrylate random copolymers containing 98.1% 
of butyl acrylate and 1.9% of acrylic acid and with vary- 
ing molecular weights were synthesized by emulsion poly- 
merization by Dow Chemical Company. Two series (A 
and B) were produced by changing the conditions of syn- 
thesis in order to obtain different architectures and, thus, 
a wide range of viscoelastic properties. Chain Trans- 
fer Agent (CTA) was used to control the weight average 
molecular weight, Mw, of the samples for a given series. 
For each sample we characterized its Mw and gel content. 
The average particle size was found to be around 380 
nm, see Table |TJ It was determined by The Dow Chemi- 
cal Company that all polymers have a low proportion of 
short branches. 
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Polymer 


CTA 


Mw 


PDI 


d 


Gel content 




(%) 


(kg/mol) 


(-) 


H 


(%) 


BglllO 




1115 


3.39 


368 


30 


A1570 




1572 


2.57 


400 




A650 


0.1 


651 


2.18 


400 





TABLE I: Properties of the model acrylic polymers. The 
parameter do is the diameter of the particles, see the main 
text for additional details. 

Latex solutions have been dried using two different 
techniques for mechanical and adhesion tests described 
below. Rheology and tensile tests required thick films 
(^600 /im), so that latexes were cast in silicone molds 
and dried during a week at room temperature followed 
by 5 min at 110 °C in an oven. For adhesion tests, thin 
films (~140 /im) coated on glass slides were made. In 
this case, latexes were cast on glass slides and dried for 
24 hours at room temperature followed by 2 min at 110 
°C in an oven. In both cases transparent cohesive films 
were obtained showing a good coalescence of the particles 
of the latex. 




FIG. 1: Processed top-view frame of the cavitation process. 
Blue contours represent the borders of the cavities our algo- 
rithm is able to detect. Only a few small cavities are missed 
because they are below the noise level. They will be tracked 
in the next frames when their area exceeds (see the main 
text). 



B. Rheology and Probe Tack Test set-up 

The characterization of the viscoelastic properties of 
the polymers was done on an ARES rheometer (TA) with 
a standard plate geometry at the Universite Catholique 
dc Louvain (UCL). The frequency range was between 
10 -2 and 10 2 rad s , while the temperature ranged from 
30 to 90 °C. We also performed tensile tests to obtain the 
mechanical properties of the sample in uniaxial deforma- 
tion. Experiments were carried out in a standard tensile 
Instron equipment (5565) equipped with a videoexten- 
someter (SVE). We imposed two different cross- head ve- 
locities v, 1.05 and 0.105 mm s _1 , for samples with an 
initial length Iq of 15 mm (initial cross section So = 2.5 
mm 2 ), resulting in a nominal initial strain rate v/Iq of 
0.07 and 0.007 Hz, respectively. 

A home built "probe tack" set up [2l| was employed 
to observe the deformation structure of the soft adhesives 
and to measure force and displacement during debond- 
ing. The apparatus consists of a cylindrical flat ended 
probe brought into contact with an adhesive layer. Af- 
ter a contact time of 10 s, the probe was pulled away at 
a constant rate V of 1 or 10 /im s" 1 . As the thickness 
of the sample h$ is 140 /im, the nominal strain rate ap- 
proximated by V/ho was 0.007 and 0.07 Hz, respectively. 
The force F and displacement d were measured during 
the whole experiment. The probe was made of stainless 
steel, a material that offers a high surface energy and 
leads to good adhesion. Furthermore, in order to obtain 
a perfectly smooth and reflective surface, the probe was 
mechanically polished. In all the experiments we used a 
probe with a diameter of 6 mm. 

A microscope was coupled to this experiment in order 



to observe the debonding structure from the top. A cam- 
era (resolution of 1292 x 964 pixels) numerical recorded 
the digitalized images. Two Zeiss lenses (1.25x and 5x) 
were used in order to get low or high magnification im- 
ages, with a field of view of 7.34 x 5.48 mm and 1.92 x 1.44 
mm, respectively. Images and force-displacement data 
were synchronized with a trigger to start simultaneously 
the probe-tack experiment and the image acquisition pro- 
cess. This trigger also controlled the frequency of the ac- 
quisition of the images, setting a frame rate of 10 and 20 
fps for a velocity of 1 and 10 (im s _1 , respectively. 



C. Image analysis 

Quantitative information about the nucleation and the 
growth of cavities can be obtained by processing the dig- 
italized top-view images acquired in probe-tack experi- 
ments. We developed a simple method to analyze these 
images by only resorting to standard routines already 
available in many packages for image processing, such as 
the Image Processing Toolbox for Matlab®. An ex- 
ample of the result of this procedure is shown in Figure 
ffl 

The algorithm detects all cavities with a surface larger 
than a threshold €a = 50 pixels. Several geometrical 
quantities, such as the center of mass, the area, the equiv- 
alent diameter and the eccentricity are measured for each 
cavity. The program also assigns an index to each cavity 
and by comparing the center of mass of cavities between 
two subsequent frames, nucleation and coalescence events 
can be tracked. For additional details see Appendix [SJ 
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FIG. 2: Storage (G') and shear (G") modulus as function of 
angular frequency (a;) for the three different materials [22J. 




FIG. 3: Nominal stress versus deformation in tensile test for 
a deformation rate A z = 0.07 (solid lines) and A z = 0.007 



(dashed lines). 



III. RESULTS AND DISCUSSION 



a markedly higher stress at large strain. 



A. Mechanical Properties 



B. Adhesion properties 



The three different materials studied differ only in ar- 
chitecture and molecular weight and the molecular inter- 
actions with a substrate should thus be the same for all 
three materials. Figure [5] shows master curves at 20°C 
of G' and G" as a function of angular frequency u> [22| . 
The curves were obtained by applying time-temperature 
superposition and it can be seen that the viscoelastic 
properties of the three materials are identical at frequen- 
cies / larger than 10 Hz. However, at low frequencies 
the rheology of A650 differs from the behavior found 
for A1570 and BglllO. The elastic modulus of A650 
decreases strongly at low frequency, leading to a mate- 
rial with a pronounced viscoelastic character. A1570 and 
BglllO on the other hand can be described as soft vis- 
coelastic solids over the whole range of frequencies. No 
terminal flow was detected for any material within the 
range of frequencies investigated. 

While linear viscoelastic properties characterize time- 
dependent relaxation processes, strain-dependent behav- 
ior is characterized using large strain properties measured 
at a given strain rate. In uniaxial extension at a fixed 
crosshcad velocity, the materials show pronounced differ- 
ences, as shown in Figure[3]by the experimental curves of 
nominal stress (Jn = F/Sq versus the deformation of the 
sample A = l(t)/lo- Macroscopic flow is observed for the 
most viscoelastic material, i.e. the A650 series, while a 
slight strain hardening behavior characterizes the BglllO 
adhesive. Although A1570 and BglllO have identical lin- 
ear viscoelastic properties at frequencies above 0.01 Hz, 
the presence of a gel fraction in the BglllO series re- 
sults in a different large-strain behavior and we observe 



Probe tack tests were carried out at two probe veloc- 
ities for the three materials. For all experiments, the 
adhesive films have an initial thickness ho and arc pulled 
by a cylindrical probe of area At ■ Experiments were re- 
peated several times and the force Ft and the displace- 
ment d = h(t) — ho as a function of time were measured. 
The nominal stress is given by 



Ft 

A~t~' 



while the nominal deformation reads 

A- k 
h 



(1) 



(2) 



and represents the nominal deformation of the whole 
sample in the vertical direction. 

The experiments are filmed at low and high magnifi- 
cation to capture the dynamics of cavity nucleation and 
growth. During the displacement of the probe the vol- 
ume between the probe and the glass slide expands. As 
the adhesive is incompressible and does not slip at the in- 
terface, this increase in volume leads to a large increase 
in tensile stress inside the layer and to the nucleation of 
cavities at the interface between the probe and the ad- 
hesive [23| and to their subsequent growth. Note that 
as the volume of the cavities expands the pressure inside 
the cavities tends towards zero. 

The nominal stress-strain curves un = /(A) are shown 
on Figure H] and are discussed together with the differ- 
ent dynamics of cavity growth. At a debonding rate of 
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FIG. 4: 

1 fim s 



Nominal stress crjv for the three materials as a function of the deformation A at a pulling velocity of 10 /im s 1 (a) and 
1 (b). 



10 //m s _1 (Figure , three different shapes of stress- 
strain curves are observed for the three materials used. 
BglllO, the most elastic material, shows a sharp stress 
peak, followed by a fast decrease of ct/v • This shape is ex- 
plained by the nucleation of cavities during the increase 
in (Tat. These cavities first expand in the bulk of the 
layer but eventually coalesce at the interface with the 
substrate. This rapid coalescence leads to the fast de- 
crease in nominal stress observed and results in interfacial 
debonding. For A1570, cavities also mainly nucleate dur- 
ing the initial increase of the nominal stress. At higher 
deformation the nominal stress is found to stabilize at 
a nearly constant value, characteristic of the growth of 
cavities in the bulk and the subsequent formation of elon- 
gated walls or fibrils. At the end, the fibrils detach from 
the surface, leading to an adhesive debonding. The ex- 
periment with A650 shows a double plateau, characteris- 
tic of liquid-like materials. In this case the walls formed 
between growing cavities are too liquid-like to sustain 
the pressure difference between the low pressure cavities 
and the atmospheric pressure and pressure equilibration 
takes place before final fibril detachment ■ In this case 
cohesive failure, i.e. residues on the probe, are observed. 

At 1 /im s _1 (Figure UJa), the shape of the stress-strain 
curve of the BglllO and A650 are qualitatively identical 
except for a decrease of the overall stress during debond- 
ing. For A1570, a transition is observed towards a liquid- 
like behavior with two plateaus. 



C. Evolution of the load-bearing area 

Due to the presence of cavities, the force applied on the 
disk-shape sample is effectively only applied on a load- 
bearing cross section that becomes increasingly smaller 
as A increases. By analyzing the projected area covered 
by the cavities and subtracting it from the initial contact 
area, the average true stress applied can be calculated 



(see the following section), instead of the nominal stress 
studied in previous investigations 0, EH EH, Hit ■ 

By means of the image analysis method described in 
the previous section we can measure the total contact 
area Ao which changes with time. Note that Aq is typ- 
ically slightly smaller than At, the area of the probe, 
even in the beginning of the experiment. The small dif- 
ference between these two areas originates from the im- 
possibility, in our setup, to illuminate the whole surface 
in contact with the probe. For each frame we measure 
the total area covered by the cavities At and then deduce 
the load-bearing cross section of our disk as a function of 
time, A e = Aq — This latter quantity is simply the 
effective area of the walls between cavities. Note that, as 
the observation direction is normal to the disk, the max- 
imal diameter of each cavity is observed in the projected 
image, see the sketch in Figure [5j 

The precise measurement of the growth dynamics of 
the cavities can unfortunately not be undertaken for the 
complete force-displacement curve. Due to loss of con- 
trast and resolution we can only precisely track cavities 
until A = 3 — 5, i.e. the first part the curves shown on 
Figure |4] and all the following results will be restricted to 
this deformation range. 

The study of the evolution of the projected areas taken 




FIG. 5: Representation of the film under deformation and of 
the disk considered in the top- view analysis. 
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FIG. 6: Evolution of the ratio Ab/Ao as function of the nominal deformation A. These seven experiments were performed at a 
constant pulling velocity of 10 fim s -1 (a) and 1 /im s _1 (b). 



by the cavities Af, as a function of time and nominal 
deformation A gives interesting insights on the average 
geometry of the cavities and can be linked to the Theo- 
logical properties of the material and to the adhesion at 
the interface with the probe. 

Cavities are observed to nucleate at different locations 
for each experiment due to the fact that nucleation de- 
pends on parameters that arc difficult to control, such as 
dust or microparticles on the surface. Also the precise 
moment of nucleation can slightly vary between exper- 
iments Interestingly, however, when the probe is 
pulled at 10 /im s" 1 the function Ab/Ao (shown on Fig- 
ure^) is very reproducible for different experiments with 
the same material and is found to be similar for the three 
materials. 

On the other hand, when the probe is pulled more 
slowly (at 1 /im s _1 ) the picture is different, see Fig- 
ure [6b. Contrary to the tests at 10 /im s _1 , some scatter 
is observed for Ab/Ao for each material and Ab/Ao now 
seems to depend on the material. BglllO shows a faster 
increase in the projected cavity area, then A1570 and 
A650. This shows that at slow pulling rate, cavities grow 
with different shapes for the different materials. Cavi- 
ties growing in the more elastic material cover more sur- 
face and are thus less elongated in the tensile direction, 
whereas cavities in the softer material show a slower in- 
crease in projected areas, indicating a more elongated 
shape. This result is consistent with what was found by 
Yamaguchi et al. [3] for adhesives with different crosslink 
densities. 

The differences in the measurements between the two 
probe velocities are interesting. In fact, they show that at 
10 /im s _1 the shape of the cavity is fully determined by 
the high frequency behavior of the materials, which docs 
not vary much between the different materials. On the 
other hand, at 1 /im s _1 , differences in rheological prop- 
erties do lead to significantly different kinematics which 



will eventually lead to very different levels of dissipated 
energy. 

D. Projected shape of cavities 

During the early stages of the debonding process de- 
scribed in the previous paragraph, the shape of the pro- 
jected area of individual cavities undergoes a transition 
from a circular to a more irregular form. Initially cavities 
grow in a circular manner. As the cavities start to occupy 
more volume they begin to feel each other through elas- 
tic interactions and viscoelastic flow. These interactions 
lead to a deviation from their initial circular shape and, 
eventually, to the coalescence of cavities, further mod- 
ifying the overall shape. A simple way to quantify this 
geometrical transition is to compute the size of the differ- 
ence between the shape of the cavity and the circle with 
the same projected area placed at the center of mass of 
the cavity, see Figure [7] This absolute difference between 
areas, Ad, provides a measure of the average change in 
shape of the cavities, thus quantifying in this way how 




FIG. 7: From left to right: i) Image of the cavity, ii) Detected 
perimeter (blue solid line) and equivalent circle (red dashed 
line) placed on the center of mass of the cavity, iii) Absolute 
difference between the two areas (black region). 
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FIG. 8: Evolution of the excess area Ad/Ao as function of the load bearing area Ab/Ao at the pulling velocity of 10 /im s 1 (a) 
and 1 jj,m s~ (b). 



the material responds to an external deformation. The 
elasticity of the material acts here like a surface tension 
and restricts sharp changes in shape [25|. 

The data is best shown as a function of the relative 
area occupied by the cavities (Ab/Ao) to avoid effects due 
to differential nucleation present for the slower pulling 
velocity. For all materials and strain rates, the projected 
area of the cavities becomes markedly non-spherical as 
cavities interact or merge with each other. 

The evolution of the normalized Ad/Ag for the two ve- 
locities is shown in Figure [8] At both velocities the more 
clastic material BglllO maintains more circular cavities 
consistent with its more elastic character. This strongly 
suggests that the level of elastic energy stored in the ma- 
terial during deformation has an effect on the curvature 
of the cavities. 



E. Growth rate of individual cavities 

We estimate the growth rate of individual cavities from 
the evolution of the projected area of each cavity as a 
function of time shortly after their nucleation. Images 
of the whole probe have not enough resolution to pro- 
vide this information and we thus use high magnification 
images (5x) of the central part of the sample. 

The increase in area of a single cavity normalized by 
the area of the probe is shown as a function of time on 
figure IH1 From this figure it is clear that the growth of 
cavities does not follow a simple functional form, in agree- 
ment with previous observations [ill, EH, EH ■ Right after 
nucleation, exponential cavity growth is observed [l8j . 
but quickly after this initial stage they start to inter- 
act with the surrounding cavities and their growth slows 
down and deviates from the exponential behavior. This 
is easily explained by the fact that cavities relax the ac- 
cumulated strews in the adhesive layer very quickly after 



their nucleation, leading to a slow down of the growth. 

We aim at capturing the first stages of cavity growth, 
as differences between different materials are expected to 
be important mainly when cavities grow independently. 
Even if cavities grow exponentially right after nucleation, 
the later stages of the growth rate can be approximated 
by a square root function and a simple exponential fit 
does not permit a clean estimation of the growth rate 
a. In fact, the time variation of the area A(t) of each 
cavity reaches a maximum in a very short time and, sub- 
sequently, it decreases. A sigmoid function S can easily 
catch this behavior of A(t): 



S 



1 



-a(i-to) 



(3) 



where a is the amplitude of S (for a = 1, S — > 1 when 
t — > oo), a is the growth rate, and to is the moment of 
maximum growth. These three parameters are estimated 
from a non-linear least-squares fit of the time derivative 
of A by using the function 



dS_ 
~dt 



-a(t-t ) 



[1 



s(t-to) 



2 • 



(4) 



A typical result of the fitting procedure is shown in Figure 

El 

A box plot of the growth rate a for the three materials 
is shown in Figure [101 We have divided the cavities in 
two groups, those that have nucleated before the force 
peak during the probe-tack test (left column) and those 
nucleated after it (right column). 

First of all, one can note that for the more elastic ma- 
terials, BglllO and A1570, most of the cavities nucle- 
ate before the maximum of the stress peak is reached 
[Tt| . For the more liquid-like material however signifi- 
cant nucleation is observed even after the stress peak has 
been reached. This phenomenon can be explained by 
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FIG. 9: Example of the evolution of the area of a growing 
cavity and its time derivative (inset) as a function of time 
and example of the fit procedure used to estimate r. Points 
are experimental data from digitalized images whereas solid 
lines correspond to fits of Eq. ((3]) and Q (inset), respectively. 
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FIG. 10: Box plot of the growth rate a for the three different 
materials at the pulling velocity of 10 fim s . Cavities have 
been divided into two groups according to their nucleation 
time, before (left boxes) and after (right boxes) the force peak. 
Percentages show the proportion of cavities for each group. 
The total number of cavities were 39 for the A650, 53 for the 
A1570, and 32 for the BglllO. The box plot is characterized 
by five-numbers summaries, i.e. the smallest observation (the 
lower horizontal line), the lower quartile (lower boundary of 
the box), the median (the line inside the box), the upper quar- 
tile (upper boundary of the box), and the largest observation 
(the upper horizontal line). We have also added the mean 
of each data set (the symbol inside the box) and outliers are 
represented by stars. 



the fact that for the low modulus of the A650 material 
the compliance of the adhesive layer quickly drops be- 
low the compliance of the apparatus leading to a sudden 
transfer of energy from the apparatus to the adhesive 
layer initiating nucleation of further cavities. This ob- 
servation might thus be apparatus dependent. The most 
interesting observation is the difference in growth rate be- 
tween the different adhesives. The most elastic material, 
BglllO, and the most liquid- like material, A650, both 
show larger growth rates with a large scatter, whereas 
the growth rate of the A1570 material is found to be 
smaller and more reproducible. For the BglllO the large 
growth rate of the projected area can be explained by the 
large amount of elastic energy stored in the elastic layer, 
leading to strong cavity growth along the interface (a 
crack propagation mechanism). For A650 the resistance 
of the material is too small to prevent bulk expansion 
of cavities, also leading to rapid growth of the projected 
area. The A1570 material seems to have the optimal 
material properties and leads to a moderate growth rate. 
The large scatter in the growth rate of cavities nucleated 
before the peak, observed for A650 and BgllOO, is most 
likely due to differential nucleation at different stress lev- 
els leading to different growth rates [ll], [H|. For A1570 
the growth rate is dominated by the viscoelasticity of the 
material leading to smaller differences in the observed 
growth rates. 



F. Effective Normal Stress 

One of the most interesting results that comes from the 
detailed analysis of the kinematics of deformation is the 
analysis of the applied force. The normal component of 
the force applied to the disk Fd is the sum of two terms. 
The first one, F mi arises from the deformation of vis- 
coelastic material, whereas the other contribution, Fp, 
is due to the work done against the atmospheric pres- 
sure to increase the volume of the low-pressure cavities 
(a suction cup effect). The force Fd will be simply esti- 
mated from the nominal force Fn and the ratio between 
the area of the illuminated region and the whole probe 
area Fd = FtA /At- 

The fraction of the measured force due to the work 
against the atmospheric pressure depends on the spatial 
distribution of the cavities on the sample. Yamaguchi 
et al. used a simple model to study the dynamics of 
debonding of an axisymmetric PSA simplified to of a 
one-dimensional problem [l9l . [2p| . Their numerical in- 
vestigations showed that, after nucleation of cavities, the 
pressure field rapidly drops to zero at the position of the 
two outermost cavities, leading to a screening effect on 
other cavities inside the PSA. This result can be easily 
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FIG. 11: Convex envelope of the region occupied by cavities 
(red solid line) with area A c . 



extended to our two-dimensional arrangement of cavities 
by considering the convex envelope of the perimeters of 
the cavities. As shown on FigurcQTJ this area A c strongly 
depends on the location of cavities and can be obtained 
from the images, so that 

F P = A c {P atm - P b ), (5) 

where Pb is the pressure inside the cavities and P a tm is 
the atmospheric pressure. As Pb is of the order of magni- 
tude of the vapor pressure, P a tm ^ Pb and equation ([5|) 
reduces to flij 

F P ~A c P atm . (6) 

Although it is obvious that this crude calculation of the 
pressure field is not accurate in the nucleation region (be- 
fore and around the force peak), it gives a good approxi- 
mation after the force peak when many cavities are grow- 
ing simultaneously in size. 

From the force F m we can then calculate the effective 
tensile component of the stress applied to the material in 
the disk 



where F m = Fd — Fp. 



G. Effective deformation 

The effective stress obtained from the previous calcu- 
lations is an average value, valid for a slice where the 
cavities have their maximal diameter (i.e. where the pro- 
jected area of the walls between cavities is minimum). 
To plot a true stress versus strain curve, we should also 
consider the local average deformation along the tensile 
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X 

FIG. 12: Local averaged deformation (A) versus global de- 
formation A for the three materials at a pulling velocity of 
10 /xm s~ . The black line is a guide for the eye with slope 
one. 



direction in the wall for a position in this slice. This can 
be obtained from the load-bearing area in a way analo- 
gous to the effective stress. The deformation A = h(t)/ho 
can be rewritten as A = R^/R(t) 2 = A$/Ab considering 
volume conservation of a stretched layer without cavity 
formation. Analogous to the correction of the nominal 
stress we now use the load bearing area to write the ef- 
fective deformation as 



<A> = (8) 

In Figure fTS] this effective average value of (A) is plot- 
ted as a function of the nominal A. The results show that 
the local deformation always exceeds the nominal one, 
suggesting a localization of the deformation in the obser- 
vation plane analogous to a necking process. The necking 
process appears to be unstable (i.e. the slope of (A) vs 
A increases with increasing A) for BglllO (crack propa- 
gation at the interface due to the stress concentration at 
the crack tip) and A650 (no strain hardening and cohe- 
sive failure) and stable for the A1570 which has the best 
PSA properties. This Figure shows well how the elon- 
gational properties of the adhesives should be optimized. 
If too much elastic energy is stored during elongation, 
stresses at the edge of the cavities cannot relax and the 
cracks coalesce at a relatively low value of A. If too little 
elastic energy is stored, the debonding geometry leads to 
necking and cohesive failure. This optimized set of prop- 
erties is consistent with the PSA design rules proposed 
by Dcplace et sd. [|| and is also in agreement with the 
observations made on the growth rates from figure ITD1 
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H. Effective stress vs effective strain curves 

We can now discuss effective stress versus effective 
strain curves as presented in Figure 1131 The initial peak 
present in the nominal stress is not observed anymore 
for the BglllO material and is much less pronounced for 
the two other materials. At 1 and 10 /mi s _1 the effec- 
tive stress for the A650 keeps decreasing after the peak 
and leads, eventually, to cohesive failure. For the in- 
termediate molecular weight (A1570) the effective stress 
decreases first and then slightly increases while the most 
interesting behavior occurs for the BglllO where the ef- 
fective stress never decreases after the peak force. One 
would expect the true stress to be much more directly 
related to the material properties and it is clear by qual- 
itatively comparing Figure [13] for example with Figure 
02 that the elasticity influences greatly how the effective 
stress varies with extension. The increase in effective 
stress for the BglllO is clearly related to the cavities ex- 
panding laterally as cracks and this increase in effective 
stress reflects the presence of a stress concentration at 
the cavity edge which leads to eventual coalescence of 
adjacent cavities and debonding. The moderate increase 
in true stress of the other two materials is characteristic 
of the extension of the walls between cavities. 

To go even further, we can finally compare the ef- 
fective stress cr e as function of (A) with the true stress 
cry = F/A(t) (which, due to incompressibility, can be cal- 
culated by <jt = Act 7v ) obtained from the tensile test (Fig- 
urejSJ. Our correction of the stress and strain values from 
the debonding experiments using the load bearing area 
is a first attempt to obtain effective stress strain curves 
that can reasonably be compared to results from mate- 
rial characterization obtained by traction experiments. 
The results of this comparison are shown on Figure [T3"k 
and!13b. Obviously the two stresses arc very different at 
values of A close to 1, since the degree of confinement is 



very high 0, [H[. However, as the deformation of the ad- 
hesive layer increases the effective stress should become 
closer to the tensile stress in uniaxial extension since the 
walls between cavities are not confined anymore. This 
is qualitatively observed in Figure 113a and 113b but one 
should keep in mind that the stress remains highly het- 
erogeneous in the foam structure and is far from being 
uniaxial. Note also that for the slow pulling speed the 
contribution of Fp is more important compared to the 
faster pulling velocity and small errors made by our ap- 
proximations might thus be more important for this case. 
The most striking difference is between the A1570 and 
BglllO where an apparently small difference in uniaxial 
constitutive behavior at this strain rate results in a com- 
pletely different debonding mechanism with a completely 
different distribution of local stress. 



IV. CONCLUSIONS 

We have carried a systematic investigation of the kine- 
matics of deformation of model thin adhesive layers 
made from acrylic pressure-sensitive-adhesives, as they 
are debonded from a flat-ended cylindrical probe at two 
different probe velocities. 

The rheological properties of the three adhesives were 
characterized in the linear viscoelastic regime and in uni- 
axial deformation until rupture at two different strain 
rates. The three adhesives were chosen to show differ- 
ences in mechanical behavior at low frequency in small 
strain and at large strain due to variable levels of molec- 
ular weight and chain branching. 

The debonding of the layer from the probe occurred 
through the nucleation and growth of cavities which then 
led to an elongated foam structure. However, the rela- 
tionship between the applied force and the nominal de- 
formation were markedly different for the three adhesives 
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representative of behaviors spanning from too liquid-like 
to too solid-like. 

The kinematics of the deformation of the layer was 
characterized by image analysis as a function of time and 
the three materials were systematically compared. The 
average shape of the cavities nucleating during debond- 
ing and the total projected area of the cavities in the 
plane of the adhesive film were characterized quantita- 
tively for all three materials at two different velocities. 
Very few differences in the overall projected area were 
observed at V = 10 ^m s _1 . However, cavities were 
more spherical projected area for the more elastic adhe- 
sive at 1 fxm s _1 while cavities were the most irregularly 
shaped for the lower molecular weight adhesive. Further- 
more an estimate of the local tensile strain in the plane 
of observation showed that the local tensile strain sys- 
tematically exceeds the nominal strain and diverges for 
the lowest molecular weight (leading to cohesive debond- 
ing) and the most elastic adhesive (leading to interfacial 
failure by crack propagation) and was only stable for the 
intermediate adhesive showing the best PSA properties. 

The kinematic information was used to calculate for 
the first time to our knowledge the effective stress as a 
function of time in the stage where cavities grow mostly 
in the plane of the film and are not yet very elongated 
in the tensile direction. While this effective stress drops 
after the peak force for the two uncrosslinkcd materials, 
it keeps increasing after the peak force for the BglllO. 
Such a qualitative difference leads to an entirely different 
debonding mechanism, with stable fibrils for the two un- 
crosslinked materials and crack coalescence for the more 
elastic BglllO. 

These results show that small differences in rheological 
properties in small and in particular large strain, lead to 
significant changes in the kinematics of deformation un- 
der the same applied boundary conditions, which then 
has a great influence on the work on debonding. This 
coupling between rheological properties and kinematics 
is a great challenge for modeling soft materials and we 
hope that our results will be the base of comparison with 
simulations of computational fluid mechanics using real- 
istic material properties. 
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Appendix A: Appendix: Image processing 

We developed an image analysis tool based on sim- 
ple thresholding conversion. The algorithm starts with a 
calibration routine before the nucleation of cavities. In 
this first step, through a trial and error procedure, we 
estimate the critical level To (with < tq < 1) for the 
conversion from greyscale to binary image. Besides, we 
determine the region of the image within which we run 
our detection routine for the cavities. This region is es- 
tablished at the beginning of the image recognition pro- 
cedure and does not evolve with time. The algorithm 
assigns to each cavity an index and manages dynami- 
cally the events of nucleation and coalescence [for details 
see point 4.]. Obviously, an empty list is created at the 
start of the procedure. Then, the algorithm repeats the 
following steps for each frame: 

1. The image is filtered with a low pass filter in order 
to reduce its noise content, we typically use a simple 
averaging over windows of size 3x3 pixels. 

2. The format of the image is converted from greyscale 
to binary according to To, that is, all the pixels with 
luminance smaller than t arc mapped to 1 (white) 
while the others to (black). 

3. All the connected regions with RICH smaller that an 
threshold ea are removed. This step is easily im- 
plemented by morphologically opening the binary 
image. 

4. The boundaries between black and white regions 
are traced and labeled with an index. The children 
of each parent object are discarded in order to avoid 
the wrong detection of small cavities inside a large 
encompassing cavity. These white spots are cre- 
ated by the unscattered light that passes through 
the cavity and is reflected back from the steel sub- 
strate. Although their position and their extension 
is related to the contact region of the cavity with 
the steel substrate, these quantities are very sen- 
sitive to many irreproducible factors, such as the 
intensity of the light, the magnification factor and 
the sample alignment. For this reason these white 
spots are not taken into account in the analysis of 
the images from whole probe experiments. 

5. For each cavity several geometrical quantities are 
measured, e.g. its center of mass, area, equivalent 
diameter, and eccentricity. 

6. By comparing the center of mass of cavities in the 
current and the previous frame, the index of each 
cavity is changed according to the list of indexes 
of the previous frame. In this step the processes of 
nucleation and coalescence of cavities are handled. 
For each new nucleated cavity a new entry in the 
list is created with a new index ut + 1, where tit 
is the largest index of the list. However, when the 
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coalescence of two or more cavities occurs, the new 
data of the coalesced cavity are assigned to the low- 
est index in the list while the entries of the other 
cavities are deleted. In this manner we are able to 



track the evolution of each cavity and record all the 
coalescence events. 
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